Zero-temperature phases for chiral magnets in three dimensions 
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We calculate the energies of various multiple-spiral spin configurations for the three-dimensional 
model of chiral magnets. The ground-state phase diagram is obtained in the space of anisotropy, 
magnetic field, and interaction parameters. Regimes with multiple-spiral, or spin crystal ground 
states of two- and three-dimensional nature are identified. The Skyrmion number responsible for 
anomalous Hall effect is examined for each spin crystal configuration. A geometric interpretation 
is given for the three-dimensional spin crystal lattice structures as the periodic arrangement of 
hedgehog and anti-hedgehog singularities. 
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I. INTRODUCTION 

Skyrmionic spin texture is a natural excitation of 
a ferromagnet in two spatial dimensions (2D) carry- 
ing nontrivial topological numbers 1,2 . Although their 
creation energy remains finite (of order of several J, 
J=ferromagnetic exchange energy) and can be excited 
thermally, an important work of Polaykov 3,4 demon- 
strated that non-topological spin-wave excitations al- 
ready destroy order at all finite temperatures, pre- 
emptying the chance of a phase transition driven by topo- 
logical defects. This is in contrast to the XY spin model 
in the same dimension, where the transition to disordered 
phase indeed occurs by proliferation of topological defects 
known as vortices 5 . 

Recent advances in the experimental investigation of 
chiral magnets 6-11 brought new insights into the role and 
importance of Skyrmions in the class of magnets char- 
acterized by both Heisenberg (J) and Dzyaloshinskii- 
Moriya (D) interactions. A most interesting finding 
of the recent experiments is that a rather small mag- 
netic field of order D 2 / J leads to the proliferation 
of Skyrmions and their crystallization into a columnar 
(3D) 6,12 ' 13 or a pancake-like (2D) 9 triangular array, de- 
pending on the dimensionality. The stability of the tri- 
angular Skyrmion lattice was successfully demonstrated 
based on Ginzburg-Landau theory 6 and on microscopic 
model calculations 14 ' 15 . 

Much less understood aspect of Skyrmion physics in 
chiral magnets is whether it is possible to realize other 
forms of Skyrme crystalline order. Indeed it was argued 
that a square array of Skyrmions and anti-Skyrmions 
would be possible energetically 14,16 ' 17 . In three dimen- 
sions, study of different 3D crystalline structures was pi- 
oneered by Binz and collaborators 18 . In this paper, we 
carry out a comprehensive study of the energetic stability 
of many different kinds of spin crystal configurations in 
three spatial dimensions, under the influence of a num- 
ber of spin anisotropy interactions and magnetic field. 
Our calculation focuses on zero-temperature energetics 
and differs from other considerations where temperature 



control also played a vital part 6 ' 16,17 . It also differs from 
Binz et al. ls in that the role of magnetic field in stabiliz- 
ing Skyrmion order is more thoroughly examined here. 

Some conclusions reached in our paper are as follows, 
(i) Two-dimensional multiple-spiral structures are ener- 
getically unfavorable in three spatial dimensions, (ii) 
For strong interactions fully three-dimensional spin crys- 
tal phase becomes more stable compared to single-spiral 
phases. Typically, magnetic field tends to stabilize 3D 
spin crystal structures, (iii) All 3D multiple-spiral struc- 
tures can be understood as a lattice of hedgehogs and 
anti-hedgehogs. Such identification allows a complicated, 
but still straightforward evaluation of the Skyrmion num- 
ber. Organization of the paper begins by introducing 
several multiple spiral states in Sec. II along with some 
single spiral states. The free energy for each state is 
worked out following the detailed procedure described in 
the Appendix. In Sec. Ill the ground state phase dia- 
gram is worked out in search for the energetic stability of 
multiple spiral phases in three dimensions. We also carry 
out a calculation of the spin chirality (Skyrmion number) 
for each of the multiple spiral phases in Sec. IV. A use- 
ful visual interpretation is given of the three-dimensional 
Skyrmion crystals as a lattice of hedgehog singularities. 
Finally we conclude with a summary in Sec. V. 

Throughout the paper, the expression "multiple spiral" 
state is used interchangeably with "spin crystal" state. A 
single spiral state does not form a spin crystal. 



II. MULTIPLE SPIRAL STATE MODEL 

The multiple spiral spin configuration can be written 
in the general form 



n + *r, 

N 



n*e- jk ° 



(1) 



as the sum of a uniform magnetization no = (0, 0, no) due 
to the Zeeman field and N independent spirals of wave 



2 



vectors k G 



1 • • • N. It is assumed that all k^'s are of Minimizing it further in terms of no gives 



the same length, and the amplitudes n a ■ n* = n 2 H /2 are 
the same for all the spiral modes. The magnitude of the 
spin is position-dependent, and we impose the average 
normalization on n r 



(n2) = l = n 2 + 2^n a -n* =n 2 +Nn 2 H . (2) 

a 

In this sense, our spins are "soft" , with the origin of 
softness possibly traced to the coupling of spins to the 
conduction electron bath or the averaging over fast fluc- 
tuations. 

For right-handed spirals we can choose 



n Q = — e °(e lQ - ie 2a ), 



(3) 



using a pair of orthonormal vectors forming a triad, e\ a x 
62a = k a . An arbitrary phase 9 a is assigned for each 
spiral mode for generality. 

Our consideration of the energy of multiple spiral spins 
is based on the free energy 



B 



n 



D 2 /J' 



E 



1 B 2 



2D 2 /J' 



(7) 



It is interesting that all multiple spiral states of the type 
(1), regardless of the number of spiral components TV 
and their orientations, carry exactly the same k-vector 
length, magnetization no, and the ground state energy, 
in the presence of magnetic field. 

It is then up to the anisotropy and interaction terms 
C, U, W to separate out the energies of different multiple 
spiral states. Below we work with eight exemplary spin 
configurations, drawn from both recent experiments and 
other theories, with A, the number of spiral components, 
ranging from 1 to 6. For A > 1 one will obtain spin 
crystal state of some sort. 

First, three kinds of N = 1 single spiral states are 
considered. 

• Helical spin state with k = (fc/v / 3)(l, 1, 1) has the 
spin configuration 



F [n] = j d?vf[n], 



f[n] = -A(n 2 x + n 2 y + n 2 z ) + Dn- (V xn)-B-n 
+ l -J[{Vn x ) 2 + (Vn y ) 2 + (Vn z ) 2 ] 



(On x y fon v y /on z y 
\dx) \ dy J \dz ) 

+ U(n 2 x +n 2 y +n 2 z ) 2 + W(ni + n y + ni). (4) 



Ferromagnetic exchange and Dzyaloshinskii-Moriya 
(DM) interactions are given by J and D while C, U, 
and W are the spin anisotropy and interaction terms up 
to fourth order introduced by Bak and Jensen 19 . Zee- 
man field B is always directed along the z-direction in 
this paper. The A-tcrm can be dropped in practice be- 
cause we will always require average magnetization of 
unity, which also makes our calculation effectively zero- 
temperature. The integration J d 3 r extends over three- 
dimensional space and is divided by the space volume. 

If we remove all but J, D, and B from the free energy, 
the multiple spiral state defined in Eq. (1) can be readily 
shown to have the energy 



E = ^Nn 2 H ( Jk 2 - 2L»fc) - Bn . (5) 
Upon optimizing k at k = D/ J it becomes 



E=- l -^Nn 2 H -Bn = l -^ T {n 2 -l)-Bn . (6) 



n^ 111 = n + nf 11 V k - r + (n^ 111 )*e- 4k - r , 
n H 111 _ n H g f x + y-2z x-y 



(8) 



This state appears as the ground state at low mag- 
netic field (B < B c ) for bulk MnSi crystal 6 " 8 . 

• Another helical state has k = {k/\/2){\, 1, 0) and 
the spin configuration 



n^ 110 = n + nf 11 V k ' r + (nf 110 )*e 



Hiio„«k-r i /„Hnn \* — ik-j 



n H t0 / „ . .y- x 

— — e z+i 



V2 



(9) 



This state appears as the low-field ground state in 
the 2D geometry 9 - 14 - 15 . 

• Conical state has k = fc(0,0, 1) and the spin con- 
figuration 



n r - no + iij e + ) e 



Co ^k-r i /„Co\*^ — ik-r 



n 



Co _ n H ab 



e l0 {x-iy). 



(10) 



It includes the ferromagnetic (FM) state as a spe- 
cial case. This is the ground state for 3D MnSi at 
fields B > B c 6 - S . 

With N — 2 and JV = 3we consider two kinds of two- 
dimensional Skyrme crystal states extended along the z- 
dircction in a columnar manner. 
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• Square-lattice Skyrmion crystal (SkX2) (N = 2) 
has the propagation vectors ki = fc(l,0, 0), k 2 = 
k(0, 1,0), and the spin configuration 



(a) SkX2 



SkX3 



n; 



n 



SkX2 
r 

SkX2 
1 

SkX2 



n 



SkX2 e ik -r 



SkX2^* e -ik Q 



— e^(y-iz), 
2 v ; 



(11) 



The spin configuration is that of half-Skyrmions 
and half-anti-Skyrmions alternating on a square 
lattice. Such a state was previously considered in 
Refs. 14,16,17. Its existence in nature remains elu- 
sive so far. 




FIG. 1: The k- vectors of five multiple-spiral spin states con- 
sidered in this paper, (a) and (b) are the k-vectors of SkX2 
and SkX3 in two dimensions, (c), (d) and (e) are the k-vectors 
of SC, FCC and BCC spin crystal states, respectively. 



Triangular-lattice SkX state (N = 3) has 
ki = fc(l,0,0), k 2 = fc(-l/2,V3/2,0), k 3 = 
k(— 1/2, — -v/3/2,0), and the spin configuration 



„FCC _ v , 

n r — n 



n 



n FCC e ik a .r 



+ (nS cc )*e-*--') ) 



n: 



SkX3 



n 



SkX3 



n: 



SkX3 



„ , \"~Y_SkX3-ik a T , („SkX3\* c -ik a iA 
n +2^V n " e + > e J' 

a 

n -f^ (z + iy), 
2 I 2 2 ' 



SkX3 _ 
*3 — 



— e . * + ,_*-,- 



(12) 



This state was found to exist in the A-phase of 
MnSi and related B20-compounds 6 ~ 8 as well as in 
the 2D thin-film sample at low temperature in the 
presence of magnetic field 9 . 

Another possible spin crystal arrangement with N = 3 
is the simple-cubic (SC) multiple-spiral spin structure 18 
which consists of ki = fc(l, 0, 0), k 2 = fc(0, 1, 0) and k 3 = 
fc(0,0,l). 



,sc 



n 



^(„f e *- + („f)*e-^ r ), 



n 



sc 



,sc 



n 



SC 



^(y-iz), 
^(z-ix), 
-^-e lBz {x - iy). 



(13) 



At N = 4 one obtains face-centered-cubic (FCC) 
multiple spiral spin structure 18 consisting of ki = 
(fcA/3)(l,l,l), k 2 = (fc/V3)(-l,-l,l), k 3 = 
(fc/\/3)(-l, 1, -1), k 4 = (fc/v / 3)(l, -1, -1), and the spin 
configuration 



fcc _ n H je, ( x + y-2z x-y 

ni ~— e { V6 

fcc _ !^je 2 ( x + y + 2z x-y 



n 9 



V2 



n. 



fcc _ n H ie 3 fx + 2y + z z-x 



+ i- 



V2 



n 



FCC _ n H 



2x + y + z y — z 
+ 1 



V2 



(14) 



Finally, at TV = 6 one obtains body-centered- 
cubic (BCC) multiple spiral states 18 consisting of 
^ = (fc/x/2)(l,l,0), k 2 = (fc/\/2)(l,-l,0), 
k 3 = (fc/x/2)(l,0,l), k 4 - (fc/\/2)(-l,0,l), k 5 = 
(fc/\/2)(0, 1, 1), k 6 = (k/V2)(0, 1,-1), and the spin con- 
figuration 



n 



BCC 



n 



n 



BCC _ n H J0 1 



BCC _ n H A0 2 



n 



n Bcc = M e « 3 



,BCC _ n H A6 



z + i 
z + i 
y + i 
y + i 



y-x 
V2 
.y+x 

V2 

x—z 

x+z 
72 



= — — e x+i 



,bcc _ n ^ t8 5 (~ , ; £-J/ 
V2 
z + y 



,BCC _ n ff 



e ° x+i 



(15) 



The last three multiple-spin configurations (SC, FCC, 
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BCC) were first considered by Binz et al. ls under the 
general effort to understand the "partial-order" phase of 
bulk MnSi at high pressure 20 . Their method of analyz- 
ing the interaction effects is different from ours, and the 
magnetic field effect was not extensively discussed in their 
work. 

By feeding all eight spin configurations into the free 
energy form (4) it becomes possible to evaluate their en- 
ergies as a function of B. The quadratic theory consisting 
of J, D, and B alone were already shown to give degener- 



ate energy, Eq. (7). Evaluation of the various anisotropy 
energies can be found in the Appendix and here we only 
quote the final expression. The final energies are conve- 
niently expressed in the reduced units, 



j^ c >Dyj^ u >DVJ^ w >Dyj^ B - (16) 



H[1H] : + u(-\nt + \nl + l) + W (-\n* + n 2 + 1) - Br*, 

H[110] : Ittcjz + u (r 2n ° + 2nl + + w (- 2 i n ° + y n ° + Bna > 

Conical : i(n 2 - 1) + U + W^-n^ - + ^)~ Bn , 

SkX2 : \{nl 1) + + §„» + 5) +w(-\n i + \nl + \) - Bn , 

mvn 1 n 2 Q -l / „ n-nl\ 3 / 2 19 4 7 2 17\ 
„ / „ /l-ng\3/2 35 4 19 2 29 \ „ 

SC : l(n 2 - 1) + u(-n 4 + ^n 2 + + \nl + ^ - Bn , 

FCC : iif^l + ^(-^ - I, 2 + g) + ^(-InS + Jn 2 + |) - Bn , 

BCC : ^ + Kr* - l»s + I) + + ^ + 1) 

/I — n 2 \ 3 / 2 

- 4V27C/ 2 + 30f7VT + 9W 2 ra ( — -Bn . (17) 



The normalization n 2 , + Nn 2 H = 1 is used to express the 
energies solely in terms of no. The above energies can 
be minimized for no between and 1 for a given field 
strength B 7 and such energies can be compared among 
the eight candidate states. In particular we are inter- 
ested in the parameter regime (C, U, W) in which one of 
the multiple-spiral states becomes the ground state, with- 
out invoking the finite-temperature effects as in previous 
Ginzburg-Landau calculations 6, 16 ' 17 . 



III. GROUND STATE PHASES 

Previous two-dimensional lattice model study yielded 
a phase diagram with the ground states H[110], SkX2, 
and SkX3 appearing over various sectors of the (C, B)- 
plane 14 (The parameter A 2 in Ref. 14 corresponds to C 
in the present model) . Our energy calculation in the con- 
tinuum model reproduces qualitatively the same phase 
diagram (Fig. 2) when the search for minimum energy 
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Fig. 3, we conclude that the two-dimensional Skyrme 
crystal phases are stable only when the corresponding 
spatial dimensionality is also two. As shown in Fig. 3, 
FCC state enjoys a particular energetic advantage in the 
high-field region while the low-field side is dominated by 
spiral states. A larger U value tends to push the FCC 
region down to a lower magnetic field [Fig. 3 (c) and (d)]. 



IV. SKYRMIONS AND HEDGEHOGS 



FIG. 2: Ground state phase diagram for two-dimensional 
spin states, {H[110], SkX2, SkX3}. In a previous work 14 
only C > was considered. Interaction parameters are 
(U, W) — (0,0.1) for this phase diagram. SP refers to full 
spin polarization, or ferromagnetic state. 



was confined to three candidate states of two-dimensional 
character, {H[110], SkX2, SkX3}. 
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FIG. 3: Ground state phase diagram for all the states we 
considered, (a) (U,W) = (0.1,0.1), (b) (U, W) = (0.1,0.5), 
(a) (U, W) = (0.5,0.1), (a) (U, W) = (0.5,0.5). 

When three single spiral states are compared in en- 
ergy, either H[lll] or conical states win out, leaving no 
room for H[110] ground state for any of the parameter 
ranges we searched. We conclude that H[110] found in 
the thin-film chiral magnet 9 was a special feature of two 
dimensionality. 

Now all the states considered in this paper are com- 
pared together and the resulting phase diagram is shown 
in Fig. 3. SkX2 and SkX3 states have disappeared, re- 
placed by spiral states H[lll] or conical at low field. Be- 
sides the (U, W) values shown, we have searched for a 
wider set of parameters but still did not find a region with 
SkX2, SkX3, or SC ground state. Comparing Fig. 2 and 



In the previous section we have established the phase 
diagram of the soft-spin model in the variational space of 
single and multiple spiral phases. The main interest of 
this paper, we emphasize, is in the identification of three- 
dimensional phases carrying a non-trivial Skyrmion num- 
ber. A general consideration relates the Skyrmion num- 
ber to the effective magnetic field felt by the underlying 
conduction electrons and the emergence of anomalous 
Hall effect (sometimes called topological Hall effect) 21 . 
With a three-dimensional space one can define such num- 
ber for each xy plane along the z-axis, the direction of 
the applied field, as 



Q{z) 



i 

47T 



, 9 fdn dn\ 
d r n • I — -x — I 



V dx dy J 



(18) 



The integration extends over one unit cell in real space 
in the xy plane. This quantity is an integer 2 for the unit- 
modulus field satisfying (n r ) 2 = 1. Our soft-spin config- 
urations generally do not obey this condition locally, but 
only on average, (n 2 ) = 1. Let's therefore improve upon 
the above definition by using the normalized vector field 
n = n/|n| to define the Skyrmion number with proper 
topological characteristics as 



Q[z) 



l 

4-7T 



1 

4-7T 



,2 - f9h 

drn id~x 
i 



d 2 r 



Oh 

dy 
fdn 

\dx 



dn 

dy 



(19) 



The second line follows from the first by a simple manip- 
ulation. In the following, we will be largely concerned 
with the evaluation of Q(z) in each of the multiple-spiral 
phases and how it depends on the applied field through 
the uniform magnetization Uq. 



A. SkX2 & SkX3 

For the two-dimensional spin crystal SkX2 the spin 
configuration, up to an overall constant, is given by 



n oc (sin y, cos x, cos y + sin x + in). 



(20) 



We choose k = 1 and m = Uq/uh for convenience. Actu- 
ally one finds a periodic array of "nodes" (localized points 
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where the magnetization vector n vanishes) for to = 0. 
Expansion around a particular such node (xo,yo) gives 
the linearized spin structure 



n oc (ycosyo, —xsinxo,m — ycosy + xsmx ). (21) 

The coordinates are measured from the nodal position. 
It leads to the Skyrmion density 



1 



P(m) 



4tt [X 2 + y 2 + P(m) 2 ] 3 / 2 ' 



where 



P(m) 



m sin xq cos yo 
\J\ - cos 2 x sin 2 y 



(22) 



(23) 



and {X, Y) is a linear coordinate transformation of {x, y). 
The node positions (xo,j/o) are ( — 7r /2,0) and (ir/2,n) 
in a unit cell of size (2tt) x (2tt), and both points give 
P(m) = —m. Integrating Eq. (22) over (X, Y) gives the 
Skyrmion number 



Q(m) 



--sgn[m]. 



(24) 



As to passes from positive to negative, the Skyrmion 
number changes from —1/2 to +1/2. With two such 
nodes per unit cell, the total Skyrmion number per unit 
cell depends on to as — sgn(m). It is convenient to view 
the situation in the enlarged three-dimensional space con- 
sisting of two spatial coordinates (x, y) and one param- 
eter to. The set of coordinates (x,y,m) = (x o ,y o ,0), 
where (x , yo) is the node position of SkX2 at m = 0, de- 
fine the centers of "hedgehogs" . As is well known, a pas- 
sage through the hedgehog center changes the Skyrmion 
number by one unit. 

One can carry out a similar analysis for SkX3, for 
which the spin structure reads 



TABLE I: Nodes of the magnetization vector in a unit cell for 
SC. P(z,m), defined in Eq. (28), associated with each zero 
can tell whether a given node is the center of a hedgehog or 
an an anti-hedgehog (17 = sin -1 (m/\/2)). 



xo 



i-v 



yo 



-f -1 
-i-'n 
- 3 -i-v 
'f + v 



+ v 



j+r) 
-j+V 



P(z,m) 



V2 - m 2 (l - m 2 )z/2 
72 - m 2 (m 2 - l)z/2 
V2-m 2 (l -m 2 )z/2 
V2-m' 2 (m 2 - l)z/2 
V2 - m 2 (l - rn 2 )z/2 
V2 - m 2 (rr 2 -l)z/2 



V2 - m 2 (l - m 2 )z/2 
V2-m 2 (m 2 - l)z/2 



One easily realizes that a node in the magnetization 
is a rare event in the two-dimensional spin structure. 
Each component of the spin vector taken equal to zero 
define a two-dimensional surface embedded in a three- 
dimensional space. For 2D spin structures such as SkX2 
and SkX3, such a surface is extended vertically along the 
z-axis and traces out a curve in the xy plane. The inter- 
action of three curves, from each component of the spin 
vector taken to zero, defines the node. The three curves 
in a plane do not cross at a single point except through 
fine tuning, such as the tuning of uq through zero in the 
SkX2 case. With both SkX2 and SkX3, the nodeless spin 
structures also happen to be topologically nontrivial and 
carry a finite Skyrmion number. 

The same reasoning applied to a generic three- 
dimensional spin crystal structure states that nodes are 
unavoidable because three independent surfaces would 
cross at isolated points. Below we show that each such 
node is the center of a hedgehog or an anti-hedgehog, a 
passage through which the Skyrmion number changes by 

+!(-!)■ 



n oc (n x ,n y ,n z ) 



K x . V3y 
n x = V3 cos -sin— — , 

x V%y\ x 

2cos-+cos— j sin-, 

n z — cos x + 2 cos — cos -). m . 

2 2 



'v 



--( 



(25) 



There are no nodes in this magnetization profile, and by 
direct evaluation one finds Q(z) = —1 for all uniform 
magnetization — 1 < m < 1. In a sense, SkX3 is a unique 
phase carrying a nontrivial Skyrmion number without 
possessing (or even proximate to) real-space singulari- 
ties in the spin structure. Indeed, recent experiments 9 ' 11 
single out the SkX3 phase as the only topologically non- 
trivial magnetic phase so far discovered in chiral magnets. 



B. SC, FCC & BCC 

Taking k = 1, the spin configuration for SC becomes 

n oc (siny + cosz,sinz+cosa;,sina;+cosj/+TO). (26) 

The three phase factors 9 a are all assumed to be zero. 
There are eight independent zeroes in a unit cell and 
their coordinates are listed in Table I. Taylor-expanding 
around a given node (xo,yo,Zo) to first order in the de- 
viation results in 



noc (y cos yo ~ z sin z o, zcos zq — x sin xq , x cos x — ysmy ). 



(27) 



7 



The coordinates (x,y, z) are measured from their zeroes. 
The triple product n • d x ii x d y n — z(cos x cos y cos z — 
sin sin yo sin zo) follows immediately, and the normal- 
ized Skyrmion density becomes 



the node positions proved difficult for general phase an- 
gles, but we do achieve some simplification with special 
choices, e.g. Q\ = 82 = and 9$ = —64 = — tt/3. In this 
case nodes occurs at 



1 P{z,m) 
4^[X 2 +Y 2 + P(z,m) 2 } 3 / 2 

where 

. z(cosa;o cos yo cos zo — sin^o sinj/o sin^o) 

P(z,m) = : =^ , 

V 1 — cos^ xo sin uq 

(28) 

and (X, Y) is again related to (x, y) through linear coor- 
dinate change. Integral over (X, Y) yields the Skyrmion 
number 

Q(z) = ^sgn[P(z,m)], (29) 

with P(z,m) for each node listed in Table I. We de- 
fine such node with Q(z = + ) — Q(z = 0~) = +1 
as the hedgehog center, and the other one with Q(z = 
+ ) — Q(z = 0~) = —1 as the anti-hedgehog center. A 
snapshot of a pair of hedgehog/anti-hedgehog configura- 
tions obtained from Eq. (27) is shown in Fig. 4. 

As can be seen from Table I, hedgehogs and anti- 
hedgehogs occur in alternating fashion along the z-axis. 
Average over a unit cell in the z-direction yields the aver- 
age Skyrmion number Q = J Q n dzQ(z)/2-K which equals 
(ji = sin _1 (TO/\/2) 



cosz\cos{x+y) — V3sin(or— y)j = 0, 

sinz^\/3cos(a; + y) — sin(a; — y)\ = 0, 
2s'mzsm(x + y) + 2coszcos(x — y) + m = 0. (31) 

One possible set of zeroes come from choosing cos z = 0, 
V3 cos(ir. + y) = sm(x-y) and 2 sin z sin(x + y) + m — 
0. There are eight possible solutions. Another pos- 
sible set is sinz = 0, cos(x + y) — -\/3 sin(a; — y) and 
2 cos z cos(x— y)+m = 0. This gives another eight zeroes. 
A third set arises from having cos(x+y) = 0, sm(x-y) = 

and sinz^\/3cos(a;+?/)— sin(x— y)j = 0. There are 16 in- 
dependent nodes that belong to the third set as listed in 
Table II. From the Table one can read off the z-integrated 
Skyrmion number. The result is in fact twice that given 
Eq. (30) for SC if one replaces r\ = sin _1 (m/v / 2) with 
t = sm - 1 (m/2 v / 2) now. The first two sets of nodes con- 
tribute zero Skyrmion number as each z-plane contains 
an equal number of hedgehogs and anti- hedgehogs. For 
BCC spin crystal phase Binz and Vishwanath obtained 
Q = -sgn(m)(2/ 7 i-)cos- 1 (|m|/2V2) for \m\ < 2^2 and 
zero otherwise 22 . Together with our calculation, a com- 
plete information is available for the magnetization (mag- 
netic field dependence) of the Skyrmion number for bulk 
spin crystal phases. 



Q 



l-iv (0 < ry < f ) 

Ufa-§) (f <»?<§) 

i] = sin -1 (m/%/2) . 
When m exceeds v2, Q(z) is identically zero. 



(a) SkX2 



(b) SkX3 



(c)SCSFCC 



(30) 




FIG. 4: Spin configuration around (a) hedgehog and (b) anti- 
hedgehog center in the SC spin crystal. Magnetization vector 
n vanishes at the center of each singularity. 

The FCC spin configuration can be analyzed with the 
similar strategy after taking k = \/3. Keeping track of 




FIG. 5: 2-integrated Skyrmion number Q for (a) SkX2, (b) 
SkX3, and (c) SC and FCC spin crystals. The FCC result is 
multiplied by 1/2 in the figure. Magnetization m is normal- 
ized to unity at the maximum value with nonzero Q. 



C. Nature of 3D Spin Crystals 

Based on the analysis of the multiple-spiral states both 
of 2D and 3D characters, we are led to suspect that re- 
alizing a topologically non-trivial (in the sense of having 
nonzero Q{z j) spin structure free of singularities is ex- 
tremely difficult in three dimensions, whereas it should 
be quite easy in two dimensions. A genuine 3D Skyrmion 
of course is without a singularity 23 , and only a lattice of 
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TABLE II: Nodes of the magnetization vector in a unit cell 
for FCC (£ = sin" 1 {m/2-^2))). 



x 


Vo 


«0 


P(z,m) 


7r/4 


5tt/4 


7tt/4 + £ 


V8- 


m 2 (4 — 


m 2 )z 


5tt/4 


7T/4 


7tt/4 + £ 


V8- 


m 2 (4 - 


ra 2 )2 


tt/4 


7T/4 


7tt/4 - £ 


V8- 


m 2 (m 2 


-4)2 


5tt/4 


57r/4 


7tt/4 - £ 


V8- 


m 2 (m 2 


-4)z 


3tt/4 


7tt/4 


5tt/4 + £ 


V8- 


m 2 (4 — 


m 2 )z 


7tt/4 


3tt/4 


5tt/4 + £ 


V8- 


m 2 (4 — 


m 2 )z 


37r/4 


3tt/4 


5tt/4 - £ 


V8- 


m 2 (m 2 


-4)2 


7tt/4 


7tt/4 


5tt/4 - £ 


V8- 


m 2 (m 2 


-4)2 


7T/4 


7r/4 


3tt/4 + £ 


V8- 


m 2 (4 — 


m 2 )2 


5tt/4 


57t/4 


3tt/4 + £ 


V8- 


m 2 (4 — 


m )2 


■w 1 A 


^TT 1 A 


•57T/4 t; 


V8- 


m 2 (m 2 


-4)2 


5tt/4 


tt/4 


3tt/4-£ 


V8- 


m 2 (m 2 


-4)2 


3tt/4 


3tt/4 


tt/4 + C 


V8- 


m 2 (4- 


m 2 )2 


7tt/4 


7tt/4 


tt/4 + C 


V8- 


m 2 (4- 


m 2 )2 


3tt/4 


7tt/4 


tt/4-£ 


V8- 


m 2 (m 2 


-4)2 


7tt/4 


3tt/4 


tt/4-C 


V8- 


m 2 (m 2 


-4)2 



such structures would qualify as a "3D Skymre crystal" . 
All the 3D states discussed in this paper and elsewhere 18 
can be instead described as a "hedgehog crystal" . 

In the hard-spin limit where n 2 = 1 is rigorously en- 
forced at the local level, we will find that a periodic array 
of singularities becomes energetically forbidding and give 
way to single-spiral structures without nodes. It is pos- 
sible, however, that the spins remain soft in real chiral 
magnets due to interaction with the metallic host that 
can "absorb" some of the spin from the local moments 16 . 
In such a case, the singularities inherent in the hedgehog 
lattice can be realized without an excessive cost in energy. 
The other exciting possibility is that the spins remain 
hard even against the metallic background, and yet some- 
how form a topologically non-trivial 3D Skyrme crys- 
tal phase such as envisioned by Klebanov and others 24 . 
Whether such a state can be energetically stable com- 
pared to hedgehog or single-spiral phases is a subject yet 
to be explored in the context of chiral magnetism. Creat- 
ing a singularity-free Skyrmion lattice structure requires 
the inclusion of higher-order harmonics (in the sense of 
Fourier decomposition of the spin structure), which also 
partly explains why our first-harmonic construction only 
yields hedgehog lattices. 



V. SUMMARY 

The ground state energies of several single- and 
multiple-spiral spin configurations are considered, on the 
basis of the Ginzburg-Landau free energy suggested by 



Bak and Jensen 19 . Phase diagrams in the space of 
anisotropy and magnetic field (C, B) are worked out for 
a variety of interaction parameters (U, W). Previous 
two-dimensional phase diagram 14 is recovered and three- 
dimensional phase diagrams are worked out anew. Sim- 
ilar efforts were made earlier by Binz et al. to carve 
out the 3D phase diagram, focusing on the stability of 
3D spin textures such as SC, FCC, and BCC, under zero 
magnetic field and using the interaction terms which are 
non-local in space. Leonov et al. 17 considered the phase 
diagram spanned by temperature and field, with two- 
dimensional spin textures SkX2 and SkX3 as candidate 
states. Our search is confined to zero-temperature, but 
covers both 2D and 3D spin textures. 

As shown in Fig. 3, magnetic field plays a crucial 
role in stabilizing Skyrmionic order over single-spiral 
states in both 2D and 3D geometries. Generically, two- 
dimensional Skyrme crystals (SkX2 and SkX3) arc fa- 
vored only if the corresponding effective space dimension- 
ality is two. For three dimensions, 3D Skyrmion struc- 
tures such as FCC and BCC are more favorable energet- 
ically. Our calculation suggests the possibility for a spin 
crystal phase of FCC symmetry very close to zero mag- 
netic field for a sufficiently strong interaction parameter 
U. 

The Skyrmion number embodied by each multiple spi- 
ral phase was examined. An array of hedgehog and anti- 
hedgehog centers were identified in the case of three- 
dimensional spin crystals. Passage through each such 
center change the Skyrmion number by ±1. Generally 
there exists a non-zero, z-integrated Skyrmion number 
for finite uniform magnetization. We showed that SkX3 is 
the only spin crystal structure without singularities and 
yet exhibiting a nonzero Skyrmion number. Constructing 
a singularity-free spin crystal state in three-dimensional 
chiral magnet is an interesting challenge for the future. 
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Appendix A: Calculation of C-anisotropy 

To calculate the C-anisotropy one should first obtain 
dn x /dx: 

^ = E(^<e^ r - ? ^«)* e - ik ° r ), 

a 

/(^) 2 - 2^(fc S) 2 <«r. (Ai) 



9 



With similar expressions from J (dn y /dy) 2 and 
/ (dn z /dz) 2 , C-anisotropy energy reads 

cE(( fc S) 2 <«)*+(^) 2 <«)*+(^) 2 <«)j- 

a 

(A2) 

The eight states have the C-anisotropy energies 



Helical(H m ) 
Hclical(Hno) 
Conical 
SkX2 
SkX3 

sc 

FCC 
BCC 



(l/6)Cn 2 H k 2 

(l/8)Cn 2 H k 2 





(3/16)Cn 2 H k 2 


(2/3)Cn 2 H k 2 
(3/8)Cn 2 H k 2 



(A3) 



Including the C-anisotropy gives the total energy for the 
various spin states 



Appendix B: Calculation of L?-energy 

Dividing up the spin configuration n = n + <fr as in 
Eq. (1), n 4 is expanded as 

n 4 = (n 2 +2n -* + * 2 ) 2 

= n 4 + 4(n .*) 2 +4n2(n .*) 
+4* 2 (n • *) + 2n^* 2 + * 4 . 



H[I1I] state gives J(n • *) 2 = J* 2 
J <& 4 = n%, and the [/-energy 

H[lll]: U^nt + ^nlnl + n^. 

H[110] state gives / (n • *) 2 = nf,n 2 H /2, / * 2 
/ <J? 4 = n A Hl and the [/-energy 

H[110] : U (n^ + 4n 2 n 2 H + n%) . 



(Bl) 



(B2) 



'Hi 



(B3) 



Conical state gives / (n • 3>) 2 = 0, / * 2 = n 2 H , J 3> 4 
riff, and the [/-energy 



Conical : U (hq + 2n 2 ,n% 



(B4) 



SkX2 state with J* (n • *) 2 = n\n 2 H , J $ 2 = 2n 2 H , and 
J $ 4 = 5n% gives the [/-energy 



Hin 
Huo 
Co 
SkX2 
SkX3 
SC . 
FCC : 
BCC 



X -n\ (jk 2 - 2/Jfc) + ^Cn 2 H k 2 

l -n 2 H ( Jfc 2 - 2L>fc) + icn^/c 2 

^n 2 H (jk 2 - 2/Jfc) 

n|(jfc 2 -2/Jfc) 

\n\(jk 2 - 2Dk) + ^-Cn 2 H k 2 



Jk 2 - 2Dk 



2n 2 H {jk 2 - 2Dk} + -Cn 2 H k 2 
3n 2 H ( Jfc 2 - 2/Jfc) + ^Cn 2 H k 2 (A4) / * 2 ("o ■ *) 



SkX2: [/ (n 4 + 8n 2 n| + 5n^) . (B5) 
There are four types of terms contributing to j $ 4 . 

• EcK • n a)( n « • K) = 2n H 

• Ei# m ( n i- n D( n ™- n m)= 2 < 

• Ez# m ( n * ' n m )(n* • <J = (1/2)^ 

• £ i# „>rO(nf • n m ) = (l/2)n% 

SkX3 state gives / (n -*) 2 = (3/2)njX* and / * 2 (n - 
<&) is non- vanishing, 



As the remaining U and W interactions do not contain 
the fc-dependence, one can minimize k for the above en- 
ergies and obtain the optimal k- vector length for each 
spin state 



( n ° ' n ki)( n k 2 ' n k 3 ) 

ki44c2-Hc 3 =0 

9n n^cos(6>i+fl 2 + fl 3 )- ( B 6) 



There are 24 combinations of the type 



/ 



(n • m)(n 2 • n 3 ) = ± no n%e i & +0 * +9 *) (B7) 



fc Hl11 

fc SkX3 
fc FCC 



D/(J+C/S), fc Hll ° = £>/(J+C/4), 

fc Co = D/J, fc SkX2 = D/J, 

D/(J+C/8), fc sc = D/J, 

D/(J+C/3), fc BCC = D/(J+C/8). (A5) 



and another 24 from their complex conjugates contribut- 
ing to / 3> 2 (n • 4>). J* 4 gets contributions from the 
following four terms. 

• EaK-<)(na-<)=3n| f 
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• Ei# m ( n * • n **)( n m • <J = 6n% 

• £i*>r °m)(n? • <J = (27/8K 

• E^ ro (nrn^)(nf • n m ) = (3/8K 
Overall, SkX3 has the {/-energy 

SkX3: C/^ + 9non^cos(6»i+0 2 + 3 ) 

+12n^+^n 4 ff ). (B8) 

The phase angles influence the total energy as <~ cos(#i+ 
02+63)- It can be minimized by choosing 6*1+^2+^3 = 7r. 

SC state gives f(n • *) 2 = n\n\ and J * 2 = 3n#. 
/ $ 4 = 12n|j is obtained from 

• E a ( n « • n a)( n <* - n a) = 3n tf 

• E^ m ( n * ■ n V( n m ■ O = 6n% 

• Ej^fo • n m )(n ; * • <J = (3/2)n% 

• E i#m (n ; -n; i )K-n m )^(3/2)n 4 ff 
{/-energy for SC reads 

SC: U (n 4 + 10n 2 n^ + 12n 4 H ) . (B9) 

FCC state gives J* (n • *) 2 = 4n 2 ) nf / /3, and / * 2 = 
4n^. / * 4 = (68/3)^ derives from 

• Ea( n «-<)( n «- n a) = H 

• Ei #m ( n * ' n D( n m ' <J = 12< 

• Ez^ m ( n * ' n m)K • <„) = (8/3)n^ 

• Ez #m ( n * ' <JK • n m ) = 4r4 

A careful calculation shows that contributions from all 
four wave vectors different give zero contribution to J <& 4 . 
The U -energy reads 

FCC : U^nU yn 2 n 2 ff + fnj^ . (BIO) 

BCC state gives / (n • *) 2 = 2n 2 Q n 2 H and / * 2 = 6n^. 
J <fr 4 = 54n|j consists of the following terms. 

• E>a-<)(n a - n*) = 6r4 

• E i#m (n ; -nr)(n m .< l )=30n 4 , 

• Ez# m ( n * ' n m )(n* • <J = 9r4 

• E i#m (n; 'OK -n m ) =9n% 



As in SkX3, 

/" * 2 (n ■ *) = J2 ( n ° • n qi )K 2 • n q 3 ) (Bll) 
qi+q2+q;?=o 

is non- vanishing from combination of vectors ki+kzt— ks = 
0, k 3 +k 6 -k 1= 0, k 2 +k 5 -k 3 = 0, k 2 +k 4 + k 6 = 0. The 
[/-energy becomes 

BCC: J/(n 4 + 14n 2 n^ + 54<) 
+ Un n 3 H ( 

5 cos(6»i +0 4 - 65) ~ V2 sin(0i + 6> 4 - f 5 ) 
+5 cos(0 3 +0 6 -0 1 )- y /2am(O s +6 6 -6 1 ) 
- 5 cos(0 2 + 65 ~ 3 ) + V2 sin(0 2 + 5 - 9 3 ) 
-5cos(0 2 + 6» 4 + 6i 6 )-V2sin(6i 2 + 6l4 + 6l 6 )). (B12) 

The phase angle dependence of the [/-energy is more 
complicated than in the case of SkX3 because there are 
six such angles instead of three. 

Appendix C: Calculation of W-energy 

Spin components are decomposed as n x = <& x , n y = <f> y 
and n z = n + & z . Hence, n 4 = $ 4 , n 4 = <I> 4 , and 

< = (n 2 + 2n 2 $ 2 + $ 2 ) 2 

= n 4 + 6n 2 ,$ 2 +4ng$ z +4n ^ + $ 4 . (CI) 

For all the states under consideration we have J n 4 = 

K 

H[lll] state gives J n 4 = (l/6)n|f and J n 4 = n 4 + 
2n^r? H + (l/6)n|f for the IF-energy 

H[lll]: W(n 4 +2n 2 n| + in 4 ,). (C2) 

H[110] state gives / n 4 = (3/32)nJj and / n\ = n 4 + 
3n^n^ + (3/8)nJj for the T^-energy 

H[110]: w(n 4 +3nln 2 H + ^n 4 H ). (C3) 

Conical state gives J n% = (3/8)n|j and J n 4 — n\ for 
the VF-energy 

Conical: wQn^ + n^j. (C4) 

SkX2 state has two types of terms in calculating J $ 4 = 
(9/4)r4. 
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• £JK)*I 4 = (3/4K 

•Y.^ m \{n l ) z \ 2 \{n m ) z \ 2 = {Z/2)n% 

With Jni = (3/8)r4 and / = n 4 , + 6n 2 ,n 2 H + {9/4:)n%, 
the W-energy becomes 

SkX2 : W(no+6nln 2 H +3n 4 H ). (C5) 
SkX3 state gives 

4n y $^ = 6no(n kl ) z (n k2 ) z (n k3 )2 

= 6n nf r cos(ei+fl 2 +«3)- (C6) 
For f$i = J^y = (81/64)n^ one works out 

•EJ(Ox| 4 = (27/64)r4 

• E^ ro IK)x| 2 |(n m )x| 2 = (27/32K 

For / $2 = (45/8)n|f one works out 

• EJK),l 4 = (9/sK 

•E^ m \(ni) z \ 2 \(n m ) z \ 2 = (9/2)n% 

The VF-energy is / n\ = (81/64)n^, J n\ = + 
6n o n%cos(6 1 +0 2 +03)+9nln 2 H + (<i5/8)n%, and 

SkX3: w(nf, + 6n n 3 H cos(e 1 +e 2 + e 3 ) 

^n 2 Q n 2 H + 2 -^n%). (C7) 

The phase angle dependence ~ cos(0i+02+03) ar i ses in 
the VF-energy as well as in the [/-energy for SkX3. 
SC state has two types of terms contributing to J $ 4 = 

• £JK)z| 4 = (3/4)r4 

• E^I(nO»l 2 l(nm)z| 2 = (3/2K 

The VF-energy is j n% = (9/4)njj, J n\ = + 6nln 2 H + 
(9/4)r4, and 



SC: w(n 4 + 6n 2 n 2 H + ^n 4 H y (C8) 

FCC state has contributions to / = f $y = J ^ 4 Z = 
(14/3)njj coming from 

• EalK)x| 4 = (2/3)n^ 

• Ez^ m \{ni) x \ 2 \{n m )x\ 2 = 

The VF-energy is / n% = (U/3)njj, J n\ = n^ + 8nln 2 H + 
(U/3)n% and 



BCC state has 



FCC: W 



(4 



•* n _ I n _ I n O 



qi+q2+q3=o 



= 3n n# 



cos(6»i + 6> 4 - 5 ) + cos(6>3 + 6 - 01 ) 



-COs(0 2 +05-03)-COs(02 + 04 + 6 )J. (CIO) 

Calculation of / $^ = j $4 = j $ 4 = (87/8)n^ can be 
obtained from 

•E Q IK)x| 4 = (9/8)n| f 

• E^ m IK),| 2 l(« m ),| 2 -(78/8K 

The VF-energy is / = (87/8)n^ and 



2 2 , 261 4 
• " II a 



BCC: w(n^ + 12nln 2 H 
+ 3n ™ff [cos (6»i + 4 - 5 ) +cos(0 3 + 6 - 01 ) 

-COs(0 2 -05-03)-COs(0 2 + 04 + 6 )]). (Cll) 

Appendix D: Minimizing the phase angles 

The three phase angles of SkX3 enter in the energy as 



{9U + 6W)n n 3 H cos(0 1 + 2 + 3 ). 



(Dl) 



It is minimized by choosing 9\ + 02 + 03 = n. 

The six phase angles appearing in the BCC energy can 
be grouped as 

X cos(0i + 04 - 5 + ai ) + cos(0 3 + 06 - 0i + a 2 ) 
+ cos(0 2 + 5 - 3 + a 3 ) + cos(0 2 + 4 + 6 + a 4 ) 

(D2) 



with X = V27t/ 2 + 30UW + 9W 2 and a t that can 
be worked out. The minimum energy configuration is 
achieved by taking all the arguments of the cosines equal 
to 7r. The energy of BCC coming from the U and W 
terms after minimizing the phase angles are 



BCC : u(n^ + Unln 2 H +5Anjj^ 

+W(n i + I2n 2 n 2 ' 261 ~ 4 



H 



8n 2 n 2 H + 14r4) . (C9) 



-4 v / 27[/ 2 + 30[/T^+9VF 2 n n^. (D3) 

No other candidate states besides SkX3 and BCC have 
the phase angles entering in the energy. 
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